Steering the potential barriers: entropic to energetic 
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We propose a new mechanism to alter the nature of the potential barriers when a biased Brownian 
particle under goes a constrained motion in narrow, periodic channel. By changing the angle of the 
external bias, the nature of the potential barriers changes from purely entropic to energetic which 
in turn effects the diffusion process in the system. At an optimum angle of the bias, the nonlinear 
mobility exhibits a striking bell-shaped behavior. Moreover, the enhancement of the scaled effective 
diffusion coefficient can be efficiently controlled by the angle of the bias. This mechanism enables the 
proper design of channel structures for transport of molecules and small particles. The approximative 
analytical predictions have been verified by precise Brownian dynamic simulations. 

PACS numbers: 05.60.Cd, 05.40.Jc, 02.50.Ey 
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I. INTRODUCTION 

The study of transport of molecules or small particles 
in narrow, confined systems such as pores, channels, or 
quasi-one-dimensional systems is of great interest due to 
it's major role in the process of catalysis, particle sepa- 
ration and dynamical characterization of these systems 
PJ-[|j]. Recent techniques enable to reveal the structural 
properties and sequence of molecules like DNA or RNA 
when they pass through narrow openings or the so-called 
bottlenecks, which cause entropic barriers [!, These 
barriers are unavoidable and present in many fundamen- 
tal processes that occur in nature, from molecules passing 
through ion channel proteins to catalytic process of ions 
in zeolites or molecular sieves. The nature of the po- 
tential barriers entirely depends on the configurational 
space available for the particle, i.e., whether the parti- 
cle moving in an one-dimensional energetic landscape or 
multi-dimensional confined geometry. In the former case, 
the nonlinear mobility increases with the noise strength 
due to noise assisted hopping over the potential 
barrier, and an opposite behavior can be observed in the 
later case [ill ]. Consequently, the transport behavior of 
the particles/molecules sensitively depends on the nature 
of the potential barrier of the system. Tuning the po- 
tential barriers is the key for controlling the transport 
of particles/molecules and allows for the development of 
new separation devices. 

Here we analyze the transport and diffusion of point 
sized particles, moving in two-dimensional narrow, con- 
fined geometries, and which are subjected to a con- 
stant bias. By changing the orientation of the bias 
the nature of the potential barriers of the system be- 
comes controllable. The dynamics of the full system 
can be approximatively described by means of an one- 



dimensional kinetic equation, the so-called Fick- Jacobs 
equation, which contains an effective potential function 
which subsumes/incorporate the geometrical restrictions 

The paper is organized as follows. In the next section 
we introduce the model system, namely the dynamics of 
a Brownian particle in a confined geometry with irregular 
boundaries. A simplified approximative one-dimensional 
kinetic description for the full system is discussed in sec- 
tion |TIIJ An analytical treatment of the problem for find- 
ing the main transport characteristics is presented in sec- 
tion [TV] The numerical techniques and main findings are 
discussed in section [V] We present our main conclusions 
in section IVT1 



II. CONFINED BROWNIAN MOTION 

The overdamped dynamics of a Brownian particle, in a 
confined geometry, subjected to the constant force F can 
be described by means of the Langevin equation written 
as 



dr 
dt 



-■F+y/lk B TZ(t) 



with 



F =Fi 



(la) 



(lb) 
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where Fii = \F\ cos(0) and F± = \F\ sin(0) denote the 
force components along and perpendicular to the 2D 
channel direction, respectively, r denotes the position of 
the particle in 2D, 7 is the friction coefficient and £(t) is 
the Gaussian white noise with zero mean and correlation 
function: (£i(t) £j (t')) = 2 <5.y S(t — t') for i, j = x,y. 

In the presence of confinement, this equation has to be 
solved by imposing reflecting (no-flow) boundary condi- 
tions at the walls of the structure. For the 2D structure 
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FIG. 1: Schematic illustration of the periodic channel confin- 
ing the movement of a Brownian particle which is subjected 
to a constant force F. L is the periodicity of the channel 
structure and 9 is the angle between the force vector and the 
positive x-axis. 



depicted in Fig. [1] the shape of the symmetric boundary 
is defined by the function 



(jj{x) — asin(27rx/ L) + b , 



(2) 



where L is the periodicity of the channel and the param- 
eters a and b control the slope and the channel width 
at the bottleneck. Here, the condition b > a should be 
satisfied in order to enable the passage of particles from 
one cell to the other. The sum and difference of the two 
parameters a + b and b — a yield half of the maximum 
width and half of the minimum width of the channel, re- 
spectively. The local width of the channel is given by 
2lu(x). For the sake of simplicity, we use dimensionless 
units. In particular, we scale lengths by the periodicity of 
the channel L, time by r = 7L 2 //cbTr, the corresponding 
characteristic diffusion time at an arbitrary but irrelevant 
reference temperature Tr , and force by Fr = jL/t. 

In dimensionless form the Langevin-equation ([TJ and 
the boundary functions © read 



dr 

dt 



= F|| e x - F x e y + VB £(t) , (3a) 
uj(x) — a sin(27ra;) + b , (3b) 

with the rescaled temperature D — T/Tr. 

III. APPROXIMATIVE ID KINETICS 

The dynamics of the system can be described by means 
of a kinetic equation, obtained originally from the 2D 
Smoluchowski equation after elimination of the y coor- 
dinate. The resulting equation, namely the Fick- Jacobs 
(FJ) equation, with a spatially dependent diffusion co- 
efficient, reads for |a/(a;)| <C 1 (the prime refers to the 
derivative with respect to x) [l5[ 



lp(x t t) = ^-V(x)e- A(x ^ D ^-e A ^ D P(x,t). 



(4) 



Here, P(x, t) is the probability distribution function 
along the length of the 2D channel, T>{x) is the spatially 



diffusion coefficient which reads 



V(x) 



D 



(1 W(x) 2 ) 



2\a 



(5) 



where a = 1/3. Therefore, the boundary function given 
by Eq. (|3b|) determines V{x). Introducing the spatially 
diffusion coefficient into the kinetic equation, improves 
the accuracy, and extends it's validity to more winding 
structures [ill, [H-[l3] • 

In Eq. 21 A(x) is the free en ergy , which in presence of 
a constant bias reads, cf. Ref. [la. Il9j. 



A(x,9) 



-Fu x- Din 



2w(i) 



sinh ( 



with 



F± u(x) 
D 



(6a) 



(6b) 



Note that the free energy A(x, 9) can be tuned by al- 
tering the aspect ratio e, which is a function of the key 
parameters, namely, the angle of the applied bias, the 
noise strength D and the shape of the periodic potential. 
Thus, it is instructive to analyze the two limiting situ- 
ations for the free energy A(x, 9) that can be obtained 
depending on the angle of the bias. 



A. Bias in the channel direction 

For 9 = 0, i.e., in the case of a constant force pointing 
in the channel direction, the dimensionless free energy 
function, Eq. [51 modifies into the form 



A(x) = -\F\x-D In [2u{x)] 



(7) 



resembling a potential function with purely entropic bax- 
riers @,[ll[231|. 

B. Bias perpendicular to the channel direction 



For 9 = 7r/2, Eq. [6]turns into 

A(x) = -D In 



2_D ginh / \F\u(x) 
\F\ m I D 



(8) 



Remarkably, this effective potential depends on the abso- 
lute value of the bias, the noise strength and the geometry 
of the structure in a non-trivial way [18( . Interestingly, 
this effective potential exhibits energetic barriers persist- 
ing even in the limit of small noise, i.e., for D — > 0. 

Thus, by changing the angle of the applied bias one can 
effectively tune the nature of the potential barriers, i.e., 
from entropic to energetic. If the direction of the applied 
bias is parallel to the direction of the channel, the free 
energy function will have a purely entropic nature, where 
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the barrier height is a function of the noise strength. At 
a finite angle of the applied bias, the free energy function 
will have both an entropic and an energetic contribution. 
At small noise strengths it is purely energetic and at high 
noise strengths it is purely entropic 



IV. TRANSPORT CHARACTERISTICS 

For calculating the key quantities of transport like the 
average particle current, or equivalently the nonlinear 
mobility, and the effective diffusion coefficient we use the 
mean first passage time (MFPT) approach. For any finite 
force the average particle current in periodic structures 
can be obtained from 

where T\(xq — > xq + 1) is the first moment of the first 
passage time distribution for reaching x — xq + 1 starting 
out at x = x . Note that the mean first passage time 
diverges for a vanishing force and consequently leads to 
a vanishing current. The effective diffusion coefficient is 
defined as the asymptotic behavior of the variance of the 
position and is given by [si— TToj] 

i2 



D eS = 



T 2 (x 



Xq 



+ l)-[T 1 (a;o->a;o + l)] 



2[Ti(x ->-x + l)]' 



(10) 



where T 2 is the second moment of the distribution of the 
first passage times which can be calculated by Eq. (fTTj) . 

The nth moment of the first passage time distribution 
can be determined recursively by means of 



T n (xo -> xo + 1) := (tn(xQ 



X 



1)) 



xa + 1 A(x,6)/D 



x T n -x(y -t x + 



d.r 



-A{y,9)/D 



1) 



dy 
(11) 



for n £ N with T n (x n -> xq + 1) = 1. Since the 
spatially dependent diffusion coefficient is periodic, i.e., 
V(x + 1) = T>(x), and A(x + 1,9) = A(x,6) - F\\, the 
recurrence relation can be simplified. The nonlinear mo- 
bility is obtained (by setting xq = 0), 



1 _ e -F l \/D 



Fu 



(12) 



I(x, 9) Ax 



and the effective diffusion coefficient is 



D. 



off 



(I 



x _-y V{x) e A ^,e)/D 



[I{z,9)YdxAz 



1 -,3 

I(x,9)dx 



(13) 



where the integral function I(x, 9) reads 

e A(x,6)/D 



Kx, 



V(x) 



A(y,8)/D 



dy. 



(14) 



A. Deterministic limit 

In the deterministic limit, i.e., D — > 0, there is no en- 
tropic contribution and the free energy function modifies 
into A(x, 9) = — Fj| x — F±cj(x), which corresponds to 
the washboard potential. The critical angle for which 
the potential barriers disappear is given by 



= arctan 



(—) 

\2ira J 



(15) 



Then, the asymptotic value of the nonlinear mobility 
reads HI HI 



lim u 



tan(9) 



tan(e c ) 







< 
> 



(16) 



and the corresponding effective diffusion coefficient reads 

1 9<9 C 

(17) 



ixm — 

D^Q D 



1 



tan(8) 



tan(8 c ) 







> 



B. Strong noise limit 

In the strong noise limit, i.e., D — > oo, the entropic 
contribution is more dominant in the free energy func- 
tion, A(x) — —F\\ x — D\n[2ui(x)}. Then, the asymptotic 
values of the nonlinear mobility and the scaled effective 
diffusion coefficient read 120, 26] 



lim /i = 



1 



D— >oo 



D f L 
———-da; / uj(y)dy 



lim 



D->oo D 



(18a) 



(18b) 



C. Zero bias limit: along the longitudinal direction 

In the absence of force component along the transversal 
direction, i.e., for the case of 9 = tt/2, one can use the 
Lifson- Jackson analytical formula for the scaled effective 
diffusion coefficient which is given by [13, [13, HH HH 



lim 



- 1 ' ^y^d, / e-^-dy 
o Jo 



(19) 



As the longitudinal force component is zero for 9 
7r/2, the nonlinear mobility vanishes. 
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FIG. 2: (Color online). The nonlinear mobility (o) and the 
effective diffusion coefficient (b) as function of the scaling pa- 
rameter D/F (F = \F\) for various angles of the force. The 
lines correspond to the analytical predictions, (a) Ea. 1121 and 
(b) Eq. 1131 whereas the symbols correspond to the results 
of the Langevin dynamic simulations for the two-dimensional 
structure with the shape defined by the dimensionless func- 
tion uj(x) = [sin(27rx) + 1.02] /2n. In limit D — > oo, both /i 
and D a ft/D reach a limiting value of 0.188 (thin black line) 
predicted by Eq. flT8j). 



V. ID VS. 2D MODELING 

Analytically derived transport characteristics have 
been compared with the results obtained from the Brown- 
ian dynamic simulations, performed by a numerical inte- 
gration of the Langevin equation, Eq. [TJ using the stan- 
dard stochastic Euler-algorithm. The average particle 
current is defined as 

(i)=lim^, (20) 
and the corresponding effective diffusion coefficient as, 

D ° S = }Too 2t ■ ^ 

For the calculations the geometrical parameters are set 
to a = 1/2-k and b = 1.02/2tt. 

The nonlinear mobility (a) and the scaled effective dif- 
fusion coefficient (b) as function of the scaling parameter 
D/F for different angles of the bias are depicted in Fig. [2] 
For 9 — 0, i.e., for a purely entropic potential, the non- 
linear mobility decreases monotonically with increasing 




FIG. 3: (Color online). Analytically obtained nonlinear mo- 
bility, Eq. (|12[) . and the scaled effective diffusion coefficient, 
Eq. (|13[) . as function of the rotation angle 6. The asymp- 
totic behaviors of fi, Eq. (|16p . and D e g/D, Eq. (|17[) . in the 
limit D — > are plotted in dashed lines. In the deterministic 
limit, both fi and D B s/D become zero for 6 > 9 C . Here, we 
present only the analytical predictions for clarity, whereas the 
agreement with numerical simulations is same as in Fig. [2] 



the scaling parameter D/F, whereas the scaled effective 
diffusion coefficient D c g exceeds the bulk diffusion D. fn- 
terestingly, for 8 — ir/4 the nonlinear mobility shows a 
bell-shaped behavior, exhibiting a peak over a range of 
D/F values, ft is due to the fact, that for very small 
D/F values the energetic nature of the barriers domi- 
nates resulting an increase in /i. Contrarily, for large 
values of the scaling parameter D / F the entropic nature 
dominates, which slows down /i. For intermediate values 
of D /F the potential function contains both the entropic 
and energetic contribution resulting a peak in fi. More- 
over, the scaled effective diffusion shows also a peculiar 
behavior for 9 = ir/4. For small D/F values, the scaled 
effective diffusion is greatly enhanced and drops down to 
the bulk diffusion constant when the scaling parameter 
D / F is increased. Note that a similar monotonic increase 
of the scaled effective diffusion coefficient with decreas- 
ing noise level was observed in periodically segmented 
channel systems [13, HI]. For 8 = tt/2, F\\ = 0, and as 
one would expect /i = 0, and D e g/D exhibits a mono- 
tonically increasing behavior with the scaling parameter, 
which could also be predicted by the Lifson-Jackson for- 
mula, Eq. (fT^|) . Note that a similar behavior of the scaled 
effective diffusion coefficient was observed in periodically 
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FIG. 4: (Color online). The numerically obtained steady- 
state distribution of particles in a two-dimensional channel 
with the shape defined by the dimensionless function io(x) = 
[sin(27ra) + f .02] /27r for different angles of the bias and at the 
scaling parameter D/F = 0.016 (F = \F\). 



segmented channel systems [23, [28| ■ 

It is worth to mention that the agreement between the 
results of the ID modeling (lines) and the numerical sim- 
ulations of the full 2D-system (symbols) is quite good, 
and the little deviation at small D/F for 9 — is a 
known failure of the F J approximation 0, H3] , where the 
assumption of equilibration along the transversal direc- 
tion does not completely hold. 

In Fig. [3] the dependence of the nonlinear mobility and 
the scaled effective diffusion coefficient on the angle of the 
applied bias is depicted. With increasing angle, the non- 
linear mobility monotonically decreases. It is the increase 
of the energetic barrier height of the potential function 
with 9 which leads to the continuous deterioration of the 
nonlinear mobility depicted in Fig. |3K- In addition, the 
diffusion coefficient shows a striking behavior: with in- 
creasing height of the energetic barriers, the scaled effec- 
tive diffusivity initially goes up, reaches a maximum at 
the critical tilt, at 9 = 9 C , then drops down. This is in ac- 
cordance to the giant enhancement of diffusion phenom- 
ena found for tilted periodic energetic potentials @, Q . 
The maximum value depend strongly on the dimension- 
less scaling parameter D/F. With decreasing D/F the 
peak of the scaled effective diffusion coefficient gets more 
pronounced, cf. the giant diffusion phenomena 0, Q . In 
the deterministic limit, i.e., for D/F — > 0, D c s/D shoots 
up at the critical angle of the rotation. But above the 
critical angle (9 C ) both the nonlinear mobility and the 
scaled effective diffusion coefficient drops down to zero, 
see Eq. QS and EI 

In other words, by tuning the angle of the bias, the 
diffusion of particles inside the channel can be enhanced 
by many orders of magnitude. This drastic enhancement 
in D e g / D was also observed in segmented channels [29[ . 



A. Steady-state distribution 

The steady-state distribution of the particles in a two- 
dimensional geometry is presented in Fig. 0] for different 
angles of the applied bias. For 9 — only the longitudinal 
force component survives, forcing the particles to move 
from one cell to the other. With increasing angle, the 
strength of the longitudinal force component decreases 
and the transversal component increases, as a result par- 
ticles are forced to move along the lower boundary of 
the channel. At higher angles (9 — > tt/2), the longitudi- 
nal force tends to zero, the particle flows in forward and 
backward direction are balanced, and as a result the net 
motion becomes zero. 
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FIG. 5: (Color online). Normalized steady-state probability 
density of particles along the propagation direction for various 
angles of the bias (F = \F\). Lines correspond to the ana- 
lytical prediction, Eq. (I22|l . and symbols correspond to the 
results of the Langevin simulations of the full system. The 
other parameters are same as in Fig. [4] 

For the reduced one-dimensional kinetic equation, 
Eq. (j4|), the normalized steady-state probability density 
P st (x,9) along the propagation direction can be derived 
analytically, and is given by [20| 



-A(x,8)/D 



P st (x, 



x+1 p A(y,9)/D 



dy 



,A{x,B)/D rx 



V(x) 



(22) 



-A(y,e)/D 



dy 



x-1 



A comparison with the numerically obtained probability 
density, is depicted in Fig. [5] for various angles of the bias. 
At 9 = only the longitudinal force component survives, 
forcing the particles move from one cell to the other, re- 
sults in a flat distribution of P s t(x, 9) over the full length 
of the cell. As the angle increases, the transversal force 
component increases and the particles are more focused 
towards the middle of the cell. For 9 = tt/2, there is no 
force component acting along the longitudinal direction, 
and particles settle down at the bottom of the cell, see 
Fig. The finite half- width in the probability distribu- 
tion is due to the presence of noise in the system. Sur- 
prisingly, the dynamics of the full system is almost cap- 



6 



tured within the approximative FJ description, whereas 
the small deviation at 9 = is a known failure of FJ 
approximation [l7l. l2(ij. 

VI. CONCLUSIONS 

We have shown that the transport nature of the Brow- 
nian particles in a channel with periodically varied cross- 
section area can effectively be tuned by altering the angle 
of the applied bias. When the direction of the applied 
bias is parallel to the direction of the channel, the effec- 
tive potential will have a purely entropic nature. In this 
situation, the nonlinear mobility decreases monotonically 
with increasing scaling parameter D/F. For 7r/4, /i 
exhibits a bell-shaped behavior which is resulting from 
the competition of the entropic and energetic contribu- 
tions to the free energy profile. For larger values of the 
angle, the energetic nature of the system rules the dy- 
namics and a monotonic increase of [i with D is observed. 



Interestingly, by rotating the bias from parallel to per- 
pendicular to the channel direction, the scaled effective 
diffusion initially increases, reaches the maximum at the 
critical angle of the bias, which depends on the ampli- 
tude of the geometric structure, and then drops down 
to zero. Summing up, by changing the angle of applied 
bias the diffusive transport of the particles can effectively 
be controlled. We have shown that the dynamics of the 
full system can be described adequately within a reduced 
one-dimensional kinetic description. 

This new approach of controlling the transport of 
Brownian particles in narrow confined structures has 
wide range of applications, including catalysis, particle 
separation, fluid mixing [311 ], and many more. 
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